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The electronic structure around an impurity in spin triplet p-wave superconductors is studied by 
the Bogoliubov-de Gennes theory on a tight-binding model, where we have chosen sinp^^+isinpy- 
wave or sin (pa; + py)+i sin (— pa;-|-py)-wave states which are considered to be candidates for the 
pairing state in Sr2Ru04. We calculate the spontaneous current and the local density of states 
around the impurity and discuss the difference between the two types of pairing. We propose that it 
is possible to discriminate the two pairing states by studying the spatial dependence of the magnetic 
field around a pair of impurities. 



PACS numbers: 74.25.Jb, 74.20.Rp, 73.20.Hb, 74.70. Pq 



The discovery of superconductivity in Sr2Ru04 by 
[ Maeno et al. jj has highly activated the field of spin- 
triplet superconductivity^- There are several key evi- 
dences suggestingthat spin-triplet pairing is realized in 
Sr2Ru04 .3, 4IElE0;ll)l3- The presence of spontaneous 
magnetic field suggested from /iSR experiment is consis- 
tent with a chiral superconductivity with broken time re- 
versal symmetry (BTRSS) 0. Tunneling spectroscopy 
' with zero bias conductance peak is also consistent 
with spin-triplet pairing with BTRSS [l^. 

Theoretically, several pairing states and/or micro- 
scopic mechanisms for spin-tri plet p-wave pairing have 

'. been proposed [H, Q Hi EE 113 • Although there are 
three bands in Sr2Ru04, it is possible to consider the 
essence of the superconducting property only by con- 
sidering quasi two-dimensional 7 -band. For simplic- 

' ity, here we classify the pairings for Sr2Ru04 into two 
types, which are qualitatively consistent with experi- 
ments of specific heat. [l8] The first type is sinpx-\-i sinpj^- 
wave proposed by Miyake and Narikiyo (MN), where 
cooper pairs are formed between nearest neighbor 
sites. The second one is sm(px-\-py)-\-isin{—px-\-py)- 
wave state proposed by Arita, Onari, Kuroki, and Aoki 
(A0KA)115|, where pairs are formed between next near- 
est neighbor sites. The pairing which has been pro- 
posed based on third order perturbation theory 0, IT^ 
has more higher harmonics where the most dominant 
component is sm(j)x+Py)+isiTi{—px+Py)-wave pairing. 
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In the following, we call sinp^j-fisinpy-wave pairing 
and sm(j)x-\-py)-\-isiii{—px-\-py)-wa.ve pairing as MN- 
type pairing and AOKA-type, respectively. The discrim- 
ination between the MN-type and the AOKA-type pair- 
ings is important because it is strongly related to the 
study of the superconducting properties and the pairing 
mechanism of Sr2Ru04. One of the remarkable differ- 
ence of the two pairing states is the winding number of 
the pairing function when we trace its phase along the 
7-band Fermi surface [l9| . The winding number is one 
in the MN-type pairing, but three in the AOKA-type 
pairing. It is interesting to propose a new idea to dis- 
criminate these two pairings having different topological 
characters. 

The aim of the present paper is to propose that exper- 
imental observation of the electronic properties and the 
internal field around the impurity can be used to discrim- 
inate the two pairing states. Actually, as shown in recent 
experiments by Lupien et al., it is possible to observe the 
local density states in Sr2Ru04 with high accuracy [20l| . 
It is shown in the previous studies that the electronic 
structure is modulated around the impurity 0, |22, 113 ■ 
For superconductors with BTRSS, circular current is in- 
duced around impurities (24] and a spontaneous field ap- 
pears by a non-magnetic impurity. In this paper, we 
calculate the local density of states (LDOS) and spon- 
taneous internal field around the non-magnetic impurity, 
and discuss the difference between the two pairings. 

To obtain the wave functions and the eigenenergies 
around the impurity, we solve the Bogoliubov-de Gennes 
(BdG) equation on a tight binding model which takes 
into account the anisotropy of the Fermi surface and the 
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superconducting gap structure [l^ . 

where Kij = —ti_j — fi6i_j and e is an index of the eigen- 
states. We set tij = t, (— 0.4t) for the transfer between 
nearest (second nearest) neighbor sites in the two dimen- 
sional square lattice of Ru atoms, and the chemical po- 
tential /i which depends on the temperature is tuned so 
that the density of electron is fixed at (n) =4/3. This 
reproduces the Fermi surface topology of the 7-sheet in 
Sr2Ru04. The energy and the temperature are scaled by 
t throughout this paper. The self-consistent condition is 

Ay = 9z,ji{{ajia^^) + (aj-^a^i)) (2) 

with (ajia^-i) = J2e'"e{^j)ueirz)f{E^) and (ajtaa) 
Y^^Ue{rj)v*{rt)f{—E^), where f{E) is the Fermi distri- 
bution function, and gz,ji is the spin triplet pairing in- 
teraction. 

The orbital part of the pair potential at each site i can 
be decomposed into sinp^,- and sinpj^-components as 

Ap^ (r,) = (A,;,,+£ - A,,,_£)/2, (3) 
Ap^ (r,) = (A,,,+5 - Au_y)/2 (4) 

as Qz^ij is non-zero only for nearest neighbor (NN) 
site pair in the MN-type pairing. For sin px±i sin py- 
wave superconductivity, we define the pair potential as 
A±{ri) = Ap^(ri)±iAp (ri). When the pair potential is 
uniform, this BdG formulation is reduced to the con- 
ventional theory for p-wave superconductors with the 
MN-type pairing functions, which has an anisotropic gap 
(sin^P2,-|-sin^Pj,)^/^. In the AOKA-type pairing case, as 
gz,ij is non-zero only for second nearest site pairs, the pair 
potential A'_j_(ri) is defined by A^^i+j^+g) and Ai^i+(_£+i;) 
instead of A^^^+i and Ai^i^y. 

To investigate the electronic structure around the im- 
purity, we calculate the LDOS at the i-th site as 

N{E, r,) = ^{|?/,(r,)p(5(£; - E,) + \v,{r,)\''S{E + E,)}.> 

We evaluate the spontaneous internal field H(r) through 
the Maxwell equation: V x H = ^j(r). The current j(r) 
is calculated as 

je{r^) = 2\e\clm{U+e,t^{a]^g,^a.i^a)} 

(7 

= 2|e|clm{fi+e,i ^[-u*(ri+e)ua(ri )/(£'„) 

a 

+ z;„(r,+e)<(r,)(I -/(£;„))]} (6) 

for the e-direction bond (e = ±x, ±y) at site r^. The 
spontaneous current je(ri) circles around the impurity 
when time reversal symmetry is broken. 

We consider a system with a square unit cell of Nr x Nr 
sites where an impurity is located at the center of the 



unit cell. By introducing the quasi-momentum of the 
Bloch state, we obtain the wave function under the pe- 
riodic boundary condition whose region covers Nk x 
unit cells. We typically consider the case Nr='i\ and 
Qz.ji = — i at a low temperature T = Q.Qlt. The chemical 
potential at the impurity site is taken as /iimp = —\Q.Qt. 
The wave function is almost zero at the impurity site. 

First, we study the properties around a single impurity, 
and compare the results between the MN-type and the 
AOKA-type pairing cases. Figure ^ shows the pair po- 
tential structure around the impurity for MN-type pair- 
ing and for AOKA-type pairing. The amplitude of the 
dominant chiral component (A+ or A^ here) vanishes at 
the impurity site, and recovers within a few site from the 
impurity, as shown in Fig. da). On the other hand, the 
other component with opposite chirality (A_ or A'_) is 
induced around the impurity as shown in Fig. H^b). The 
phase structure of the induced component, shown in Fig. 
^c), has a remarkable difference between the two pair- 
ing cases, reflecting the winding along the Fermi surface. 
In the MN case, ArgA_ has -|-2-winding at the impu- 
rity site, and — 1-winding at four sites located eight sites 
away from the impurity along the vertical and the hor- 
izontal directions. In the AOKA case, ArgA'_ has — 6- 
winding at the impurity site, and -|-l-windings at four 
sites located on the diagonal directions in addition to 
-|-l-windings at four sites on the vertical and horizontal 
directions. Since the amplitude of the induced compo- 
nent vanishes at these winding center sites, |A^| is more 
suppressed around the impurity with —6- winding, and 
the tails of |A'|_| have a shape extending toward eight di- 
rections far from the impurity. For the difference of the 
winding structure of the induced chiral component, there 
appears some differences between the MN-type and the 
AOKA-type pairings in the structure around impurities. 

In Fig. I3a), we plot the spontaneous current je(r) 
around the impurity. The current circles around the im- 
purity, and the staggered current spreads toward the di- 
agonal directions. The current in the diagonal direction 
is larger in the MN case. We show the internal magnetic 
field in Fig. |2Ib). The spontaneous field at the impurity 
site is oriented to the —z direction. This induced field 
is larger in the AOKA case, since the the spontaneous 
current around the impurity is larger [Fig. |2ta)] due to 
the — 6-winding of A'_ at the impurity [Fig. H^c)]. To 
compensate the spontaneous field at the impurity site, 
the spontaneous field is oriented toward the ~\-z direc- 
tion near the impurity along the horizontal and vertical 
directions. 

Figure|2Ia) shows the LDOS N{E,ri) around the im- 
purity. The impurity is located at site (16,16). We plot 
the N{E,ri) at the nearest site — (15, 16), at the sec- 
ond nearest site (15,15), and at the farthest site (1,1). 
The energy spectrum of the uniform state with p-wave 
anisotropic gap is reproduced at (1,1). Splitted LDOS 
peaks appear within the gap at (15,15) and (15,16) near 
the impurity. The side peak in the positive energy range 
is larger than that in the negative range for both pair- 
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FIG. 1: (Color) Spatial structure of the pairing potential for 
the MN-type (left panels) and the AOKA-type (right panels) 
pairings. We shown a region consisting of 21 x 21 sites with the 
impurity located at the center, (a) Amplitude of the dominant 
chiral component, |A+(r)| or |AY(r)|. (b) Amplitude of the 
induced opposite chiral component, |A_(r)| or |A'_(r)|. (c) 
Phase structure of the the induced component, ArgA_(r) or 
ArgA'_(r). The phase of the dominant chiral component is 
almost uniform. 
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FIG. 2: (Color) The spontaneous current je(i') (a) and the 
internal magnetic field distribution H{r) (b) around a single 
impurity for the MN-type (left panels) and the AOKA-type 
(right panels) pairings. We show a region consisting of 11 x 11 
sites with the impurity at the center. Solid circles in (a) 
indicate the impurity site. 
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FIG. 3: (Color) (a) Spectrum of the LDOS N{E, r) at sites 
(15,15), (15,16) and (1,1). (b) Density plot of the LDOS 
N{Er^O, r) at the peak energy of the spectrum near £ ~ 
within a unit cell of 31 x 31 sites. The left and right panels are, 
respectively, for the MN-type and the AOKA-type pairings. 



ings. The peak at (15,15) is larger (smaller) than that at 
(15,16) in AOKA-type (MN-type) pairing. The LDOS at 
this peak energy is shown in Fig. Olb). The tail struc- 
ture far from the impurity is similar, but the LDOS at 
nearest and second nearest sites is different between MN- 
type and AOKA-type pairings. In the case of AOKA- 
type (MN-type) pairing between diagonal (vertical) site 
quasiparticles, the low energy state is likely to appear 
in the diagonal (vertical) site next to the impurity site. 
The anisotropic gap structure on the Fermi surface also 
contributes to the difference of the LDOS. In the AOKA- 
type pairing, the gap amplitude has local minimum at 
[110] directions on the Fermi surface in addition to the 
minimum at [100] directions, as shown in Fig. 1 of rcf. 

As we have seen in Fig. |2Ib), a drastic qualitative 
difference is not seen in the internal field between the 
two pairings in the case of single impurity. However, in 
multiple impurity case, there are some cases producing 
qualitative differences between the two pairing cases in 
the spontaneous field distribution around the impurities 
due to the interference of the phase of the induced oppo- 
site chiral component (i.e. A_(r) or A'_(r) in Fig. 
As an example of this effect, we report the case when two 
impurities are located at next nearest neighboring sites, 
i.e., sites (15,15) and (16,16) in our calculation. As in 
the single impurity case, the zero energy LDOS is large 
at the nearest site of the impurities for MN-type, and 
at the next nearest site for AOKA-type. As shown in 
Fig. ^ the spontaneous current and the field distribu- 
tion is not a simple summation of the distribution of the 
single impurity case, since the strong interference is at 
work between the two impurities. It is remarkable that 
for MN-type pairing, the field intensity is strong in the 
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FIG. 4: (Color) The spontaneous current >(r) (a) and the 
internal magnetic field H{r) (b) for the MN-type (left panels) 
and the AOKA-type (right panels) pairings when a pair of 
impurities is situated in the diagonal direction, e show a 
region of 10 x 10 sites for (a) and a region of 20 x 20 sites 
for (b) around the impurities. Solid circles in (a) indicate the 
impurity sites. 

direction perpendicular to the direction in which the im- 
purities are aligned. The field orientation is toward the 
+z direction in this case. On the other hand, in AOKA- 
type pairing, the spontaneous field, oriented toward the 
—z direction, has strong intensity in the direction paral- 



lel to the impurity alignment. This is a good example of 
phenomena in which a difference in the phase windings 
of the pairing function on the Fermi surface results in 
experimentally observable quantities. 

In summary, we have calculated, on the basis of BdG 
equation, the spontaneous field and the electronic struc- 
ture around the impurities in spin triplet chiral p-wave 
superconductor, where we compare the results between 
the MN-type pairing smpx+is'mpy and the AOKA-type 
pairing sin {p^+py)+ism {-Px+Py}- 

We have shown that differences between the two pair- 
ings appear in the LDOS structure at the sites nearest 
and next nearest to the impurity site. Namely, the zero 
energy LDOS at the nearest site is found to be larger 
(smaller) than the LDOS at the second nearest site for 
MN-type (AOKA-type) pairing. When a pair of impu- 
rities is situated in the diagonal direction, the sponta- 
neous magnetic field distribution around the impurities 
is quite different between the two pairings. The sponta- 
neous field has strong intensity along the line perpendicu- 
lar (parallel) to the direction in which the two impurities 
are aligned in the MN-type (AOKA-type) pairing. We 
expect that the properties around the impurities shown 
here provide information for distinguishing the pairing 
state of Sr2Ru04. 
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